Free Collisions in a Microgravity Many-Particle Experiment. I. 
Dust Aggregate Sticking at Low Velocities 



R. Weidling% C. Guttler^'\ J. Blum^ 

^Institut fUr Geophysik und extraterrestrische Physik, Technische Universitdt zu Braunschweig, 
Mendelssohnstr. 3, D-38106 Braunschweig, Germany 
^Department of Earth and Planetary Sciences, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan 



Abstract 



o 

<N 
O 

O 
o 



Oh 
O 



Over the past years the processes involved in the growth of planetesimals have extensively been studied in the laboratory. Based 
on these experiments, a dust-aggregate collision model was developed upon which computer simulations were based to evaluate 
how big protoplanetary dust aggregates can grow and to analyze which kinds of collisions are relevant in the solar nebula and are 
worth further studies in the laboratory. The sticking threshold velocity of millimeter-sized dust aggregates is one such critical value 
that had so far only theoretically been derived, as the relevant velocities could not be reached in the laboratory. We developed a 
microgravity experiment that allows us for the first time to study free collisions of mm-sized dust aggregates down to velocities 
of ~ 0.1 cms"' to assess this part of the protoplanetary dust evolution model. Here, we present the results of 125 free collisions 
between dust aggregates of 0.5 to 2 mm diameter. Seven collisions with velocities between 0.2 and 3 cm s"' led to sticking, 
suggesting a transition from perfect sticking to perfect bouncing with a certain sticking probability instead of a sharp velocity 
threshold. We developed a model to explain the physical processes involved in dust-aggregate sticking, derived dynamical material 
properties of the dust aggregates from the results of the collisions, and deduced the velocity below which dust aggregates always 
stick. For millimeter-sized porous dust aggregates this velocity is 8 • 10"^ m s"^ 
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The formation of planetesimals, the precursors of planets, 
is initiated by the collisional coagulation of small dust parti- 
cles and aggregates. Velocities are induced from the interac- 
tion of these dust particles with the thin gas of the protoplane- 
tary disk (PPD) and the gravitational interaction with the cen- 
tral star ( [Weidenschilling and Cuzzi 1993| ). For the initially 
micrometer- sized dust grains, the collision velocities are small 
enough to let them stick to each other and form larger, frac- 
tal dust aggregate s (|Wurm and Blum", "1998" 'Blu m et al.| |2Q0Qt 
[Blum and Wum il|2000tlKrause and Blum, 2004). However, as 
the particles are getting bigger, they decouple more efficiently 
from the surrounding gas. This leads to an increase in their 
relative collision velocities ( [Weidenschilling [ |1977| ) and it is a 
priori not clear how large particles can grow by direct sticking. 
Many laboratory experiments have shown that millimeter- sized 
particles do not stick to each other if they collide at a velocity 
which is expected to occur under the conditions in a PPD (see 
review by IBlum and Wurm| ( |2008| ) and recent experiments by 
[Beitz et al. (20111)). Indeed, the sticking velocities for those 
dust-aggregate sizes are so small that they are not yet known 
and need to be studied, and that is the goal of this paper. 



1.1. The current collision model 



Our current knowledge on the outcome of dust-aggregate col- 
lisions has been largely shaped by the collision model of Guttler] 
et al. ( ,2010j. Based o n available results from laboratory exper- 



iments, |Gtittler et al.| quantified the outcome of a collision be- 
tween two dust aggregates in terms of sticking, bouncing, and 
fragmentation. Moreover, their model predicts which of these 
outcomes actually occurs for a given set of collision param- 
eters (dust-aggregate masses, dust-aggregate porosities, colli- 
sion velocity) over a wide range of dust-aggregate masses and 
collision velocities. According to this model, the sticking ve- 
locity for millimeter-sized, porous dust aggregates (i.e. a mass 
of 0.1 mg) is as small as 10"^ ms"\ which is slower than the 
expected collision velocities for these dust aggregates in a PPD 
( [Weidenschilling and Cuzzi 199 3). The model of Guttler et al. 



was implemented into a local growth simulation (Zsom et al. 



2010[ ), and the result was that the growth in a minimum mass 
solar nebula stalls at masses of approximately 1 mg. Instead of 
further sticking to each other, the dust aggregates rebound and 



are compacted as observed in the experiments o f [Weidling et al. 
("2009). The picture changed when [Zsom et al.[(^2011) included 



the vertical dimension of the disk and sedimentation of the dust 
aggregates, because particles with diff'erent growth timescales 
get turbulently mixed and the size distribution becomes wider 
all over the disk. 

Those results are expected to be very sensitive to the exact 
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velocity is merely based on theoretical assumptions as there had 
been no experiments for millimeter- sized dust aggregates at the 
relevant velocities at that time. Therefore, it is most desirable 
to directly measure the threshold velocity for small, porous dust 
aggregates and to understand the physical processes involved. 



An implicit assumption of the Guttler et al. model is that the 
dust aggregates are homogeneous (but porous) spheres. This 
means that the model is not necessarily useable for fractal dust 
aggregates. An exact value for the fractal dimension of dust ag- 
gregates in this evolutionary phase remains highly speculative 
at this point ( [Qrmel et al!] |20Q7t |Suyama et aT) |2008| |Qkuzumi| 



et al. 2009) but it is likely that dust aggregates crossing the 
threshold from sticking to bouncing for the first time possess a 
fractal dimension less than three. A measurement of the stick- 
ing threshold for fractal dust aggregates of the same size would 
be an ambitious next step but is not within the scope of this 
paper. 

An additional challenge in this context are the diff'ering ex- 
isting definitions for fractal dust aggregates, especially for the 
radius and the volume. [Jones | ( [20TT] ) recently presented a model 
parameterizing fractal and porous particles made up of finite- 
sized constituents. He describes such particles by means of an 
inflation and a dimensionality, which are derived from unam- 
biguous properties of the particles, like the largest spatial ex- 
tent of the particle or the volume the solid matter occupies. An 
advantage of this model is that particles do not have to be frac- 
tal in the strict mathematical sense to be described with these 
parameters. 

To complement the picture on collisional grain growth in pro- 
toplanetary disks, we would like to mention that other efl'ects 
are considered to play a role but are so far only sparsely studied. 
These involve magnetic fields and magnetized dust aggregates 
dNubold and Glass meief 2000 : Domini k and Niiboldl [2002 ) as 
well as electrostatic eff'ects to enhance the cross- sectional area 
and the sticking efficiency ( [Ivlev et aL 2002|) and leading to 
reaccretion of bouncing dust grains ( Blum 2004| ). Additionally, 
gas eff'ects might lead to the reaccretion of small fragments after 
a fragmenting collision ( |Wurm et"aLl 2001a|b| ). Another eff'ect 
that increases the complexity is the rotation of the dust aggre- 
gates. This does not only lead to a more complicated treatment 
of the velocity (translational and rotational) in each individual 
collision with a possible influence on the sticking efficiency, 
but may also change the resulting shape of dust aggregates in 



the first growth phase (jPaszun and Dominik 2006 ). 



1.2. Concept and background of our experiment 

In order to investigate the transition from sticking to bounc- 
ing collisions for mm-sized porous dust aggregates, the parti- 
cles have to collide with velocities of millimeters per second. 
With currently available techniques, these velocities are not fea- 
sible under standard laboratory conditions, whereas they can be 
achieved in microgravity experiments. Previous experiments 
by HeiBelmann et al. pOlO ) describe a method to achieve very 
low collision velocities with solid particles: they injected an en- 
semble of 1 cm diameter glass beads from two opposing sides 
into a flat box under microgravity conditions. In this granular- 
gas experiment, each inelastic collision between two particles 



resulted in the dissipation of energy and, thus, in a lowering 
of the kinetic energy and velocity. Each additional collision 
slowed down the particles further, lowering the average colli- 
sion velocity over time. 

For dust particles, the coefl&cient of restitution - the ratio of 
the relative velocity of two colliding particles after and before a 
colHsion - is around s ^ 0.2 ( [Blum and Miinchj [T9931 [HeiBeP] 
mann et aL|[2007| ). Compared to the average coeflficient of resti- 
tution of the glass beads of 0.64 ( [HeiBelmann et al. 2010[ ) this 
means that an ensemble of dust particles is slowed down even 
more efficiently, providing that the average collision frequency 
is similar. This is evident from the mean velocity of an ideal 
system of equal-sized particles with velocity-independent s: 



v(t) = < — -{-(I -s)no-t 

.^0 



(1) 



(HeiBelmann et al. 2010[ ) generally referred to as Haff^'s law 
( ,Haff.,,1983i r Here, vq is the initial particle velocity at time 
^ = 0, and n and cr are the number density and the collision 
cross section of the particles. 

In order to be able to follow the trajectory of each particle 



throughout the whole experiment, HeiBelmann et al. chose the 
dimensions of the chamber containing the glass beads to be 
1.5 times as high as a single sphere diameter, preventing them 
from obscuring each other. A side eff'ect of this quasi two- 
dimensional setup is that the particles collide with the walls 
very often. In the case of glass particles and glass walls this 
eff'ect can be neglected (the coeflftcient of restitution in particle- 
wall collisions in their case was significantly higher than in 
particle-particle collisions). However, porous dust particles 
tend to stick to glass walls even at moderate velocities, which 
requires a larger test volume compared to the particle size. This 
leads to a slightly diff'erent implementation of the concept out- 
lined by [HeiBelmann et aL] and the details of our setup are de- 
scribed in the next section. 



2. Experimental Setup 



Our experimental setups are described in detail in Sects. 2.1 



2.2 and 2.3 



To account for technical advances in the past 
and future, we will refer to the experiment by the acronym 
MEDEA (Microgravity Experiment on Dust Environments in 
Astrophysics) together with a number for the respective version 
of the experiment. In future eff'orts, the improvements of the 
experimental setup will help us to address some issues more 
accurately than we are able to do in this work. 



In contrast to [HeiBelmann et al.[ we chose a cylindrical ge- 
ometry for our test chambers (see Fig. [T]). We use glass vacuum 
chambers with a diameter of 25 mm and a height of 50 mm. 
These can be agitated in a sinusoidal oscillation along the cylin- 
der axis to excite the dust particles inside. In microgravity, the 
dust aggregates are observed with high-speed cameras in back- 
light illumination. According to Equation ([T]), a high number 
density of dust particles is desirable, which however constrains 
the observability of the dust aggregates. To optimize collision 
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time and observability, we chose an optical depth of approx- 
imately 0.3, resulting in number densities from ^ = 5 • 10^ to 
5-10^ m"^ depending on the particle size. For 0.75 mm diameter 
dust aggregates this leads to a collision time r = (ncrv)~^ 0.16 
s, with V = 0.1 ms"^ being a typical relative velocity at the be- 
ginning of the experiment. 

2.7. First Setup (MEDEA-I) 

The first experimental setup, which we will refer to as 
MEDEA-I, was used in an experiment campaign at the Bremen 
drop tower in June 2010. Utilizing the catapult, we obtained 
a microgravity time of nearly 9.5 seconds in each of our 10 
flights. We used four identical setups, i.e. four test chambers 
and four cameras, in each flight, resulting in 40 experiments to 
extensively test and develop the setup and to vary parameters 
like dust-aggregate size, dust-aggregate material, number den- 
sity, agitation frequency, etc. 

As mentioned above, each vacuum chamber has a cylindrical 
shape with a diameter of 25 mm, a height of 50 mm and consists 
of glass. The bottom flange is made of aluminum and contains a 
magnetic coil able to inject a permanent magnet into the cham- 
ber. The magnet serves as a fast solid particle that can break up 
agglomerates of dust particles or accelerate a small part of the 
dust particles in the chamber when colliding with them. The 
upper flange consists of a metallic grid (mesh size 100 yum), 
which allows for the evacuation of the chamber while keeping 
the dust aggregates in the chamber. On top of each chamber, 
vacuum bellows provide the connection to the vacuum system. 
The flanges and the vacuum chamber can be removed from the 
setup in order to place the particles in the chamber and clean it 
(Fig. [T] left and top middle). 

Below the vacuum chamber, a shaking mechanism enabled 
us to excite the vacuum chamber (Fig. [T] left and bottom mid- 
dle). The mechanism consists of a DC motor driving an eccen- 
tric wheel with a horizontal axis. The bottom of the lower flange 
is pushed against the top point of the wheel, which causes it to 
move up- and downwards. Springs between the baseplate and 
the bottom flange as well as the vacuum bellows ensure the con- 
tact with the wheel even in microgravity. By this mechanism, 
the test chamber was forced to oscillate in an approximately si- 
nusoidal motion with an amplitude of A = 1 mm to A = 2 mm 
and an angular frequency between oj = 271-4 and oj = 271-16 
s"^ At the beginning of the experiment, the shaking mecha- 
nism distributes the particles in the vacuum chamber, and later 
it ensures that even in case of residual accelerations that cause 
particles to drift to the top or bottom, the particles get pushed 
back to the center of the chamber which is keeping the density 
homogeneous. According to the amplitude and frequency of 
the agitated top and bottom flange, collisions with the flanges 
occurred at a flange velocity of Aoj = 0.02. .0.2 ms"^ and an 
acceleration of Aoj^ = 0.6..20 m s"^. 

Each experiment is observed with a high-speed camera op- 
erated at 500 frames per second with a resolution of 500 x 500 
pixels. In front of the experiment, a mirror with a 90° -angle 
facing towards the camera serves as a beam splitter. A mir- 
ror at each side reflects the beam towards the vacuum chamber. 



enabling us to get two rectangular views of the experiment. Be- 
hind the chamber, two LED arrays and difl'users are placed to 
illuminate the experiment (Fig. [T] right). With this system, the 
optical resolution of the camera was 120 jjm per pixel. 

All four experiments were connected to the same vacuum 
system with a vent line out of the drop capsule. Through this 
vent line, the vacuum chambers were slowly evacuated, to- 
gether with the drop tower, down to a final pressure of 10 Pa. 
The valve was then closed right before firing the catapult and 
launching the experiment. 

2.2. Second Setup (MEDEA-II) 

The setup introduced above was used for the results pre- 
sented in this paper. However, as we also qualitatively checked 
our results with an improved setup and will present more of 
these in the future, we will also address the improvements of 
our setup after the first campaign. We already conducted a 
second drop tower campaign with a slightly modified setup in 
November 2010 - again with four experiment chambers, but 
this time with a total of 15 catapult flights. We will refer to this 
setup and experiment campaign as MEDEA-II. In addition to 
varying particle size, particle density and shaker modes, a broad 
variety of materials were tested in this campaign (see forthcom- 
ing papers). 

One important change was to modify the viewing angle of the 
camera as compared to the first campaign. As our first analysis 
of the data showed, the two views separated by 90° are impossi- 
ble to correlate in this many-particle experiment. Therefore, we 
reduced the viewing angle to 30°, putting this into practice by 
substituting the beam splitting mirror in Fig. [l]by a refraction 
prism (130° x 25° x 25°) and omitting the two outer mirrors. A 
view through the prism directly yields two images separated by 
30° and we chose to look onto the wide plane of the prism with 
the 130° corner pointing to the experiment chamber. With this 
flat viewing angle it was sufficient to use only one illumination 
array per experiment, placed in a line with camera, prism, and 
vacuum chamber. 

To ensure a better performance of the shakers at low fre- 
quencies, we installed a servo controller for each DC motor. 
In contrast to a simple power line, the controllers were able to 
compensate higher mechanical loads by a higher current. Thus, 
we could operate the shaking mechanisms at lower frequencies 
without risking them to jam. This is particularly useful if no 
residual accelerations cause the particles to drift, but one still 
wants to have a means of removing them from the top and bot- 
tom walls. 

When we used the reduced pressure in the drop tower for 
the MEDEA-I campaign (10 Pa), we found that gas drag still 
played a minor role for long particle trajectories and for small 
particles (cf. [Appendix A| ). To avoid this, we used an addi- 
tional turbomolecular pump in the MEDEA-II campaign. This 
was operated while the experiment was prepared to be shot in 
the catapult before the vacuum was sealed with a valve and the 
pump flooded. From sealing the vacuum until launching the ex- 
periment it took about 30 seconds and we achieved an experi- 
ment pressure of 0. 1 Pa. Three of our four experiment chambers 
were connected to this vacuum system while a fourth chamber 
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Figure 1 : (color online) Left: The MEDEA setup consists of a shaking mechanism with an eccentric wheel and a glass vacuum chamber, which is connected to the 
vacuum system. Middle: Photograph of the setup and a close-up of the shaking mechanism. Right: Optical system of the MEDEA-I experiment with high-speed 
camera, beam splitter, and back-light illumination. The red dashed line and green dash-dotted line illustrate the optical paths. 



was not evacuated but operated at ambient pressure (10^ Pa). 
We used this setup to study free collisions between solid parti- 
cles, where gas drag plays a negligible role. 



2.3. Third Setup (MEDEA-III) 

In a third drop tower campaign with 5 catapult flights in Au- 
gust 201 1 we mainly used the setup mentioned above with only 
slight improvements including (i) an improved version of the 
excitation mechanism that is less prone to get stuck, ( ii) a better 
illumination system, fitted to the 30° optics, and ( Hi) a restricted 
experiment volume in some experiments. Fitting the size of the 
camera images, the insets to reduce this volume had a depth 
of only a few particle diameters, allowing us to increase the 
number density of the dust aggregates while keeping the optical 
depth low enough to observe all the particles. Some of the inner 
walls of the insets were coated with a nano-material to reduce 
the eff'ective surface and, thus, to prevent sticking of the dust 
aggregates to the glass walls. We will refer to the experiments 
conducted in that campaign as MEDEA-III. 



2.4. Dust analog material 

For the experiments described in the following, we used 
dust particles consisting of irregular- shaped, poly disperse Si02, 
with ^^80% of the particles having a diameter of 1 to 5 yum 
(^99% between 0.5 and 10 yum, manufacturer infor mation) and 
a material density of 2600 kg m"^ (see Table 1 in |Blum et aL 
([2006) for more details on its properties). Dust from the stor- 
age container forms natural aggregates, which were sifted be- 
fore the experiments to obtain a narrow dust-aggregate size dis- 
tribution (see inset in Fig. [2]). We used sieves of 0.5, 1.0, and 
1.6 mm mesh sizes, which resulted in particles with sizes either 
between 0.5 and 1.5 mm (henceforth referred to as 1 mm dust 
aggregates) or between 1.0 and 2.0 mm (hereafter referred to as 
1.5 mm dust aggregates). The largest particles can always be 
slightly bigger than the mesh size due to the fact that the parti- 
cles are ellipsoids and slip through the grid if their semi-minor 
axes are smaller than the mesh size. 

To measure the volume filling factor (the fraction of the 
dust-aggregate's volume actually occupied by dust particles), 
we examined dust aggregates with diameters of about 0.5 to 4 
mm from the storage container. We developed a setup in which 
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Figure 2: The volume filling factor of the dust aggregates used in the MEDEA- 
I experiments. On average they possess a filling factor of 0.35, denoted by 
the dashed line. The gray bars denote the errors due to uncertainties in the 
measurement of the mass and volume of the particles and the dotted lines show 
the standard deviation from the mean value. The inset shows a photograph of 
typical dust particles as used in the experiments described in Sect. [3] 



a particle can be placed on a rotatable platform with a vertical 
rotation axis. Each individual dust aggregate was observed with 
a high resolution camera (1280 x 1024 pixels, 4 jim per pixel) 
and back-light illumination. After a full rotation of the particle 
around one axis with about 100 images taken, we measured the 
diameter of the particle for each horizontal line in the acquired 
image sequence. By determining the largest and smallest di- 
ameter we then calculated the volume of each slice, assuming 
an elliptic shape. Adding up these volumes and accounting for 
shadow effects of the sample holder, this yields the total volume 
of a particle. Although a polygon would be a better approxima- 
tion for each horizontal slice, this would overestimate the parti- 
cles volume if shadowing effects due to concave regions occur. 
As the particles in general have a convex shape, the assumption 
of elliptic slices, where only information about the largest and 
smallest diameter are required, is more accurate in our case. 

After the determination of the volume, the dust aggregates 
with masses ranging from 0.1 to 30 mg were weighed to an ac- 
curacy of 0.1 mg. The results are presented in Fig. [2j where a 
larger error in the filling factor is evident for the smaller parti- 
cles, which is mainly due to the uncertainty in the mass deter- 
mination. In contrast to that, the error in the determination of 
the volume is very small and can barely be seen. As the mea- 
surements for the volume filling factor scatter equally to higher 
and lower values around the average and all of the dust aggre- 
gates are formed in the same way, we assume the same volume 
filling factor of = 0.35 ± 0.05 (standard deviation) for dust 
aggregates of all sizes. The shapes of our dust aggregates can 
also be expressed by an inflation / according to the model of 



Jones] pOTT] ). As the dust aggregates consist of monomer parti- 
cles with a size distribution, which do not occupy positions on 
a regular grid, we follow his equation (34) and arrive at / = 2.9, 
which is in this case the same as the enlargement parameter de- 



fined by |Qrmel et aLl ( [20Q7] ). 

In order to determine the irregularity of the MEDEA-I parti- 
cles, we measured the aspect ratio of our particles by analyzing 
the change of the cross section of rotating particles in the mi- 
crogravity experiments. If we follow the particle trajectories 
over a sufficient number of rotations, we can assume to have 
observed the dust aggregates from all sides and, thus, also ob- 
served the largest and smallest axis. Assuming that the inter- 
mediate axis is observed in both cases, the ratio of the largest 
to the smallest cross section equals the ratio of the largest to 
the smallest axis. In our case, a typical aspect ratio is l:x:1.7, 
where 1 < x < 1.7. Using equation (33) of |Jones| ( |2011| ), this 
yields a dimensionality of 2.2 for rod-like particles to 2.6 for 
disk-like dust aggregates. 

3. Results 

We completed the analysis of one experiment with 1.0 mm 
dust particles and two experiments with 1.5 mm dust parti- 
cles from the first drop tower campaign in June 2010 with the 



MEDEA- 1 setup (cf. Sect. 2.1 ). Other experiments with similar 
sized dust particles (including the measurement campaigns de- 



scribed in Sect. |2.2| and [23] ) were so far only visually inspected 
but qualitatively agree with the results presented here. 

In order to analyze the experiments, the images from the 
high-speed camera were first separated into the two perpendic- 
ular views and then the images were binarized. The threshold 
was selected such that the size (i.e. area) of the particles in the 
center of the image was reproduced as well as possible in the 
binary image, as we used the cross section to calculate the parti- 
cle mass assuming spherical shape and the volume filling factor 



presented in Sect. |2.4| Due to imperfect illumination, particles 
at the darker edges appeared slightly larger in the binarized im- 
age, but this has no influence on the detection of the center of 
mass. 

In the next step, the position of each particle that did not over- 
lap with another particle and could be distinguished from the 
background was tracked for as long as it was visible. The posi- 
tion was computed as the center of the area in the binary image 
and then stored together with the frame number (i.e. time) and 
the cross sectional area. To achieve statistically unbiased re- 
sults, we followed every single particle in the experiment with 
the 1 mm particles consecutively. Every time a particle could 
unambiguously be identified, it was tracked as long as it could 
be followed. We then returned to the image where this particle 
was first tracked and continued to look for the next particle from 
there. For the larger 1.5 mm dust aggregates, we did not find 
any sticking event and only aimed at measuring the collision 
velocities for a representative number of bouncing collisions. 
We looked for these collisions and only followed the colliding 
particles before the incident. In the following, we will focus on 
the analysis of the 1 mm dust aggregates if not stated otherwise. 

If a particle track deviated from the ballistic trajectory, it was 
visually checked whether a collision with another dust parti- 
cle or the chamber wall had occurred. In the former case, the 
number of the tracked particle, the time of the collision and the 
duration of the undisturbed track before and after the collision 
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were noted. Collisions with the wall were ignored, as they slow 
down the particles but are of no other interest for us. To de- 
duce the particle velocity, we linearly fitted the uninterrupted 
parts of the tracks. For tracks longer than 60 sequential images, 
we split the track into segments as a track may not be perfectly 
linear due to minor gas drag eff'ects (cf. [Appendix A\. The 



velocity of each of these segments was calculated separately 
for the X- (perpendicular to the shaking direction) and for the 
z-component (parallel to the shaking direction), then geometri- 
cally added up and assigned to the point in time in the middle of 
the segment. We did not succeed in correlating the two camera 
views for maintaining full three-dimensional trajectories. This 
is why we changed the setup in the MEDEA-II and -III cam- 



paign to a 30° viewing angle presented in Sect. 2.2 in which the 



two images are in a closer relation to each other. In Appendix 



^we show that the velocity component parallel to the shaking 
direction is about twice as fast as the components perpendicu- 
lar to it. If we choose a z-axis along the direction of excitation 
(perpendicular to the top and bottom flanges), there is no reason 
why the velocities in the orthogonal x- and j-directions should 
be difl'erent to each other. This means that the dimension we 
do not see is comparable to the slower one that we do observe. 
Adding up the components geometrically this results in an un- 
derestimation of the velocity of 13 % if we only take the value 
from the projection. Nevertheless, we will use the projected ve- 
locity henceforth, as we only know a statistical but not the exact 
value for the dimension that is not observed. 

In Fig. [3] the absolute two-dimensional velocities of all 
tracked 1 mm dust aggregates during the experiment are plotted 
over time. Each dot represents one dust aggregate in one time 
interval. At the beginning, the particles move with velocities 
of 0.1 to 0.01 ms"\ which they achieve from collisions with 
the top and bottom flange. Two distinct events are marked with 
vertical lines: after t = 1.6 s, the shaker was turned ofl' and the 
mean velocity immediately decreases due to inelastic collisions 
of the particles. Although energy was also dissipated before, the 
shaking induced additional energy which kept the average ve- 
locity approximately at a constant level. Some particles are also 
decelerated in this phase, as after about 0.8 seconds the cloud of 
particles reaches the upper flange, is reflected and the number of 
collisions increases dramatically. At the second event, after 3.5 
s, a magnet was accelerated into the vacuum chamber (cf. Sect. 
2.1 ), introducing a fast, solid object colliding with several of the 
dust aggregates on its way. These particles were accelerated to 
high velocities, but mutual collisions quickly slowed down the 
dust aggregates again. Some particles appear to move at veloci- 
ties between 10"^ and 10"^ m s"\ where the lower value of this 
range is given by our resolution limit. In reality, all of these 
seemingly extremely slow particles are sticking to the wall of 
the vacuum chamber and a residual motion of the chamber or 
the particles, resulting in slight changes in the binary image of 
the particles, is then misinterpreted as a very low velocity. As 
we will in the following only treat free collisions between dust 
aggregates, these particles are not considered any further. Most 
of the free-floating particles possessed velocities between 10""^ 
and 10"^ ms"^ by the end of the experiment, resulting in colli- 
sions with relative velocities also in this range. 




Figure 3: The average velocity of the 1 mm dust aggregates decreases from 
about 0.1 ms~^ to a few millimeters per second during the experiment. The 
vertical lines denote the times where the shaker was switched off and then a 
fast solid object was injected into the experiment chamber. 



It is obvious from the velocity evolution that the dust aggre- 
gates collide with one another. However, most of these colli- 
sions could not be analyzed as particles may be colliding within 
dense clusters or only while close to other particles (overlap- 
ping in the binary images). However, a total number of 103 col- 
lisions between 1 mm dust aggregates could be observed well 
enough to analyze them in more detail. We do not expect these 
to be biased in any way and therefore regard them as a repre- 
sentative subset. For the 1.5 mm dust aggregates, we observed 
22 representative collisions. Further analysis of these collisions 
is presented below. 

3.1. Collisional outcomes 

The most prominent result of each collision is its outcome 
in terms of sticking, bouncing, and fragmentation. From the 
103 analyzed collisions of 1 mm dust aggregates, one resulted 
in the fragmentation of one of the dust aggregates, seven in co- 
agulation and 95 in bouncing. Two exemplary sequences for a 
sticking and a bouncing collision are presented in Fig. |4] After 
a sticking collision, the two particles rotate around their com- 
mon center of mass and, as a requirement for positive detection 
of sticking, we demand that this rotation is at least 180 degrees 
if we do not have any other evidence that they actually stick 
together, like information from the other viewing direction or 
a collision efl'ecting both particles (see Sect. |4.3| ). In the top 
part of Fig. |4j two dust aggregates collide with a relative ve- 
locity of 9 • 10"^ ms"\ stick together and are rotating around 
each other several times. The image sequence in the bottom of 
Fig. |4]shows two dust aggregates colliding with 6.2 • 10"^ m s"^ 
which bounce ofl' after the contact (at ^ = 100 ms), noticeable 
by the growing distance between them afterwards. At ^ = and 



^ = 50 ms, the injected magnet described in Sect. |2.1| can be 
seen in the bottom right corner. 

All of the 22 collisions of the larger 1.5 mm dust aggregates 
resulted in bouncing. Although not all particles were tracked in 
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Figure 4: (color online) Top; An image sequence of a sticking collision at a 
velocity of 9 mms~^. After sticking to each other, the particles rotate around 
their common center of mass. Bottom: For higher velocities (62 mm s~^ in this 
case) the collisions often lead to rebound. At f =100 ms, the particles are in 
contact and clearly separate on the subsequent images. The unequal choice of 
the time steps is due to shadowing effects of other particles in the many-particle 
system. A movie of this collision can be found in the online version of this 
article. 

these experiments, the analyzed colHsions are a representative 
subset, as we could not detect any collision resulting in sticking 
or fragmentation. 

In Fig. [5] all analyzed collisions are arranged according to 
their relative velocity and the mass of the smaller dust aggre- 
gate in the collision; the colors denote sticking (green), bounc- 
ing (yellow), and fragmentation (red), the symbol distinguishes 
between the experiment with the small 1 mm dust aggregates 
(circles) and the two experiments with the larger 1.5 mm par- 
ticles (squares). The color in the background is the expected 
collisional outcome following the model of Guttler et al. (2010) . 
The majority of our data matches the model as it predicts bounc- 
ing for the velocities we observed in our experiment. One par- 
ticularly fast collision, triggered by the injection of the magnet, 
occurred at a relative velocity of 1.7 m s"^ and resulted in frag- 
mentation of a particle. 

However, we also observed for the first time millimeter- sized 
dust aggregates stick to one another under controlled experi- 
mental conditions. The collision velocity in these cases ranged 




0.001 0.010 0.100 1.000 

velocity [m/s] 

Figure 5: (color online) In this velocity-mass parameter plot, most of our colli- 
sional outcomes agree with the collision types predicted by |Guttler et al.H2010] 
background color). Green color denotes sticking collisions, yellow describes 
bouncing, and red marks fragmentation; the circles represent the results for the 
1 mm dust aggregates, the squares those for the 1.5 mm dust aggregates. All of 
the collisions that led to sticking in the experiment are in disagreement with the 
model of a sharp transition from sticking to bouncing. 

from 3 • 10"^ to 2 • 10"^ ms"^ This was unexpected as the 
dust-aggregate collision model predicts sticking of 0.1 mg dust 
aggregates to occur only at velocities lower than 10"^ ms"^ 
The fastest collision leading to sticking has to be considered 
apart from the others, though. As can be seen in Fig. [5] the 
smaller dust aggregate had a mass of only 10"^ g, leading to a 
mass ratio of the collision partners of 133. In this case, a small 
fragment stuck to another particle, violating the assumption of 
the model that the particles are roughly equal sized. Therefore, 
we neglect this collision in our further considerations. 

From our experiments, we can conclude that the change from 
sticking to bouncing is not abrupt but is rather characterized 
by a transition zone where particles stick or bounce with some 
probability. This will be discussed in detail in Sect. |4j but 
before we will statistically characterize the collisions in further 
detail. 

3.2. Impact parameters 

The impact parameter h describes how centered or off'-center 
a collision is taking place. It is a measure for the vector com- 
ponent of the distance between the centers of mass of the two 
colliding particles at the impact time perpendicular to their rel- 
ative velocity vector. By normalizing it with R = ri + r2, where 
ri and r2 are the radii of the two colliding particles (assuming 
spherical shape), a value of b/R = indicates a perfectly central 
collision while b/R = I denotes a grazing collision. 

The positions of the centers of mass were calculated from the 
fitted trajectories of the 1 mm dust aggregates before and after 
each collision. In the cases where the particle tracks after the 
collision could not be analyzed (i.e. for sticking collisions or 
in collisions where one dust aggregate was partly covered), the 
existing data was insufficient to calculate the impact parameter 
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properly. Therefore, we had to discard 31 colHsions, leaving a 
total of 72 collisions including 6 sticking collisions, which we 
examined individually. Figure [6] shows the distribution of the 
normalized impact parameter bjR, marked by the crosses, in a 
cumulative diagram. The six sticking collisions are marked by 
arrows. Within the limits of the small number of events, they 
can be regarded as being randomly distributed over the impact 
parameter range, which shows that the impact parameter is not 
the main driver for the sticking in these collisions. The im- 
pact parameters of all analyzed collisions also seem to be ran- 
domly distributed (gray line in Figure [6]), but the expectation 
for the distribution of a two-dimensional projection of impact 
parameters is rather represented by the black solid line, which 
was deduced from a simple Monte Carlo simulation of arbi- 
trary collisions. The reason why central collisions are expected 
to be more dominant in the projection is that the cross section in 
the direction that is not observed is much bigger than for graz- 
ing collisions. From the comparison of our data with the black 
dashed line, representing the same distribution but with a dif- 
ferent normalization, it is obvious that large impact parameters 
near unity are over-represented. There are several possible rea- 
sons for this: (i) with the dust aggregates being elliptic rather 
than perfectly spherical and the normalization of the impact pa- 
rameter being based on an equivalent spherical radius, impact 
parameters b/R> I are possible. Although we do not see any of 
those, this still leads to an overestimation of the impact param- 
eter. ( ii) Irregular surfaces may counteract against the argument 
of decreasing cross section with increasing projected b/R (see 
above). (Hi) As shown in Fig. [7] the coefficient of restitution 
is generally larger with bigger impact parameter, meaning that 
less translational energy is lost in the collision and the particles 
separate faster than they would after a central collision. The 
less time (i.e. the fewer images) passes after a collision until 
the particles can be tracked again, the more likely it is that no 
other dust aggregate has moved to the same area in the mean- 
time and the collision partners can still be observed. 

The coefficient of restitution £ is a measure of the loss of 
energy in a collision. It is defined as the ratio of the relative 
velocity of the colliding objects after and before a collision. 
While e = denotes a perfectly inelastic collision, a value of 
£ = I means that no energy was lost in the collision, i.e. a 
perfectly elastic collision. In Fig. |7]the squared coefficient of 
restitution of 66 bouncing collisions in which we could follow 
both particles before and after the impact is plotted over the 
squared impact parameter. This analysis is limited to the 1 mm 
dust aggregates because we did not follow the rebound veloci- 
ties in the experiments with the larger particles. Although our 
2D experiment shows a wide scatter in the coefficient of resti- 
tution, a general trend of a higher coefficient of restitution with 
growing impact parameter can be observed. The linear Pearson 
correlation coefficient of the squared values is 0.5, indicating at 
least a moderate correlation. [Blum and Munch| ( p^993 ) showed 
this trend to be linear and derived a value of = 0.51 for 
the squared coefficient of restitution for a grazing collision (i.e. 
b/R = 1) of spherical particles with rough surfaces that pre- 
vent sliding during contact. Using a chi- square fit we obtain 
values ofs^ib^lR^ = 0) = 0.12 ± 0.04 and e^{\) = 0.51 ± 0.12, 
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Figure 6: The distribution of the impact parameter of the 72 investigated dust- 
aggregate collisions in the two-dimensional view. The solid gray line shows 
a random distribution of 2D impact parameters, the black solid line represents 
the expected distribution for a two-dimensional viewing geometry in a three- 
dimensional experiment, and the black dashed line is the best fitting distribution 
taking the collision statistics into account. The six collisions leading to sticking 
are marked with arrows and appear to be randomly distributed. 



in agreement with the analysis of [Blum and Mtincti| ( |1993| ), al- 
though our particles are not perfectly spherical and are likely to 
have small diff'erences in porosity and size. 

As a reason for the scatter in the coefficient of restitution 
data, we identify the rough and irregular surface and the un- 
known rotation of the dust aggregates. As an example, a dust 
aggregate can store rotational energy in one collision and re- 
lease it into translational kinetic energy in the next collision. 
Moreover, the model of |Blum and Munch] ( |1993| ) assumes 
spherical particles, which yield a linear relation between the 
impact parameter and the sine of the impact angle. In the first 
approximation, our particles are represented by ellipsoids, for 
which the diff'erence in impact angle can easily be of the order 
of 15° for any given impact parameter. Moreover, the size ratio 
of our particles can be up to 1:3 in extreme cases, while it is 
assumed to be unity in the model. 



4. Discussion 

In this section we will discuss the actual thresholds between 
sticking and bouncing collisions. 



In Sect. 4.1 we will present 
our results in terms of a velocity-dependent sticking probability 
and in Sect. |4.2| we adapt a theoretical collision model to deduce 
the velocity threshold below which all dust aggregates stick. In 
Sect. |4.3| we will discuss a peculiar collision between a single 
dust aggregate and two dust aggregates sticking together. From 
this we will get insight into the binding forces between sticking 
dust aggregates in our experiments. A brief discussion of the 
implications of these results on the growth model is given in 
Sect.l44l 
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Figure 7: The two-dimensional coefficients of restitution of 66 bouncing col- 
lisions between 1 mm dust aggregates show a significant scatter but possess 
a general increasing trend with increasing impact parameter Although only 
weakly correlated for the squared values as plotted here (linear Pearson corre- 
lation coefficient of 0.5), a linear fit agrees with the results of |Blum and Munch| 
(T993) . 



Figure 8: The transition between sticking and bouncing collisions is not sharp 
but follows a probability distribution in velocity (solid line). The black crosses 
denote the fraction of sticking collisions in the given intervals. The vertical ticks 
at and 1 denote the individual bouncing and sticking collisions, respectively. 
The asterisk marks the sticking velocity Vstkk = 0-21 mms~^ derived in Sect. 



4.1. Sticking probability 

From the results of the experiment with the smaller 1 mm 



dust aggregates presented in Sect. 3.1 it is apparent that stick- 
ing of particles also occurs at velocities higher than proposed 
by the dust-aggregate collision model of [Guttler et al.| ( [2QT0| ). 
Although we could observe only seven collisions in which the 
particles stuck to one another, which is just a small fraction of 
the total number of observed collisions, there is a clear trend 
of sticking happening at lower velocities. For the following 
considerations we discard the single sticking event of the two 
different- sized dust aggregates. However, this collision does 
not change the overall picture. 

As we observed sticking and bouncing in the same range of 
velocities and impact parameters, we propose that there is no 
sharp threshold between sticking and bouncing collisions, but 
rather a smooth transition over roughly two orders of magni- 
tude in velocity. Figure [5] shows the probability of a collision 
leading to sticking for five velocity intervals between v = 10"^ 
ms"^ and v = 0.2 ms"^ (black crosses). The probability was 
calculated by dividing the number of sticking collisions in the 
interval by the total number of collisions. Single collisions are 
indicated by the vertical bars at sticking probabilities (bounc- 
ing) and 1 (sticking). The solid line in Fig. [8] shows a logarith- 
mic fit 



/^stick(v) = -0.6 - 0.4 . Ig ■ 



ms~ 



(2) 



(for 8 • 10"^ ms"^< v < 2.6 • 10"^ ms"^) to the three values of 
the sticking probability that are nonzero. Although this proba- 
bility should be taken with care, it nicely reproduces the 100 % 
probability sticking velocity of 2.1 • 10"^ ms"\ which we will 



derive in Sect. 14.21 

In the collisions between 1.5 mm dust aggregates we did not 
observe any sticking events. To statistically estimate the num- 



ber of sticking collisions to be expected for these dust aggre- 
gates, we use the velocity of each bouncing collision to cal- 
culate the expected sticking probability. Adding up all these 
probabilities yields an expected number of 2.2 sticking colli- 
sions in these experiments, which should be observed if the 1.5 
mm dust aggregates behaved similar to the 1 mm dust aggre- 
gates. Correcting the velocity by a factor of (mi/ 1712)'^ = 1.7, 
with mi and m2 being the average mass of the smaller colli- 
sion partner of the 1 and 1.5 mm dust aggregates as displayed 
in Fig. [5] respectively, leads to an expected 1.2 sticking colli- 
sions. This correction takes a mass dependency of the sticking 
probability into account following the model by IGiittler et al. 
and Sect. |4.2| However, we do not see any sticking collision 
for the 1.5 mm dust aggregates. Using Eq. ^ we calculated 
the probability of observing no sticking collision for the given 
number of collisions with their respective velocities to be 8%. 
If we correct the velocities as above to account for the higher 
masses, we get a probability of 25% of not seeing any stick- 
ing collision. Therefore, our observations do not contradict the 
expected result. 

4.2. A model for the sticking velocity 



3.2 



From the coefficients of restitution presented in Sect 
we can learn important dynamic material parameters of the 
millimeter-sized dust aggregates. [Thornton and Ning (1998) 
developed a model for the velocity dependence of the coeflS- 
cient of restitution in collisions between elastic-plastic adhe- 
sive spheres, described by their Equations (80) and (81). In 
Fig. [9] we present our coefficients of restitution as a function of 
the impact velocity (black crosses) and the mean values of 1 1 
consecutive data points (red squares; arithmetic mean in e, ge- 
ometric mean in v). There is a lot of scatter that we attribute to 
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Figure 9: (color online) Our data for the coefficient of restitution as a function 
of collision velocity for all impact parameters. The red squares are averages 
over 1 1 single values, which are compared to the model by [Thornton and Nmg] 
jl998]f , denoted by the green solid and dashed curves. 



the distribution of impact parameters (see Figs. [6] and [7]), irreg- 
ular surfaces, and rotation. Values above e = 1 are possible if 
dust aggregates are rotating before the collision and transform 
rotational energy into translational energy; one collision even 
had a coefficient of restitution of £ = 2.1. In the averaged data, 
however, a clear trend of decreasing coefficient of restitution 
with increasing velocity can be found, which can be compared 
to the model by Thornton and Ning Here, we ignore that the 
model was developed for central collisions while our collision 
parameters are randomly distributed. We will comment on that 
point later, but if we keep this constraint in mind, we will arrive 
at interesting results. The only free parameters in the model 
are the sticking velocity Vstick and the yield velocity Vyieid of 
the dust aggregates. The yield velocity is the minimal velocity 
for which the so called yield pressure ;?yieid is achieved at the 
contact surface, which marks the onset of plastic deformation. 
The presented curves in Fig. [9] are mainly determined by the 
yield velocity, the sticking velocity defines a cut-off' at low ve- 
locities, which occurs in this plot at velocities lower than 10"^ 
ms"^ Our mean data (red squares) agrees with the model for 

+ 1.8 e-1 

-0.6 



V'yield 



0.9^n^ mm s ^ (green solid and dashed curves). The 



^ = 910 kgm"^, we can use Equation ^ to calculate Young's 
modulus of our dust aggregates to be £" = 8100 Pa, which is 
about an order of magnitude higher than earlier assumptions 
(e.g. IGiittler et al. |2009^ and oth ers). The sticking velocity in 
the model of Thornton and Ning] (their Equation (54)) is 



Vstick - 4.2 



5 4 



1/6 



(4) 



where yeif is the surface energy of both surfaces in contact, 
which we call eff'ective surface energy. This surface energy is 
a complicated parameter in this context as we have two very 



complex surfaces in contact. The model by Thornton and Ning 
assumes a smooth contact surface that uniformly contributes to 
the contact energy. In our case, we use the effective surface 
energy as a combination of the known surface energy (for a 
smooth surface), the dust aggregate porosity, and the Hertz fac- 
tor (the ratio between contact surface and cross section of two 
monomer grains) and apply 



Teif = r • ' 



^2 • 



(5) 



The first term on the rhs. of Equation ([5]) is the surface energy 
of the Si 02 material measure d for micrometer- sized monomer 
grains by Heim et al. (1999) to be y = 0.037 Jm"^ (energy 



of both surfaces in contact). The second term is the volume 
filling factor, which we use as a surface filling factor of = 
0.35 and that accounts for the dust- aggregate porosity. Still, 
two grains in contact do not have a contact which is as large as 
their cross section but smaller than this, due to the Hertz factor. 
The Hertz factor, (a/ao)^, is the ratio between the contact cross 
section of two spherical dust monomers na^ and the monomer 
cross section TOq. We approximate the contact geometry of our 
irregular grains by that of spherical grains (the geometric mean 
diameter of our grains is 1.0 yum. Young's modulus is 41 GPa, 
and Poisson's ratio is v = 0. 17) to calculate the Hertz factor. For 
the contact radius a we assume spheres sticking to each other 
only due to their adhesion. The contact radius is then given by 
Johnson et al. ( |1971 their Equation (20)) as 



a = 



97rr(l - y^)af^ 



relation between yi eld pressure and yield velocity as given by and the Hertz factor is consequently 



Thornton and Ning follows from their Equation (10) as 



v^yield 



= 6.4' 



Vyieid' 



2^ /97rr(l -v^)f^^ 



(3) 



(6) 



(7) 



where r, m, and E are the radii, masses and Young's moduli of 
similar particles, respectively. 

As a measure for the yield pressure, we use the static com- 
pression measurements by |Blum et al. ( 2006| ), who used dust 
aggregates from the same dust as in our experiments. For a 
pressure of 200 Pa, the dust-aggregate material started to yield 
to the applied pressure, which we therefore take as the yield 
pressure. If we assume dust aggregates with r = 0.5 mm and 



For the grains assumed here, this results in a^/a^ = 1.4 • 10"^ 
and the eff'ective surface energy including both geometric ef- 
fects is as low as yeif = 1-8 • 10"^ Jm"^. Figure 10 graphi- 
cally illustrates the results of the above considerations: based 
on Equation ^ we intend to deduce the sticking velocity for a 
material with sparsely known Young's modulus and an eff'ective 
surface energy which largely depends on geometry. The best 
estimate for Young's modulus is 8100 Pa (see above), which is 
the black, diagonal line. The surface energies are denoted by 
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Figure 10: (color online) Graphical illustration of the sticking velocity follow- 
ing Equation |4j for a range of surface energies and Young's moduli. Our best 
values for the surface energy (1.8 • 10"^ Jm-2) and Young's modulus (8100 Pa) 
result in a sticking velocity of 0.21 mm s~^ . 



the green vertical lines and the best value for the sticking veloc- 
ity is at 0.21 mm s"^ which is represented by the intersection of 
the left vertical line (best value for the effective surface energy) 
and the black diagonal line (best value for the dust- aggregates' 
Young's modulus). A wide range of Young's moduli from the 
softest up to solid materials (gray lines) as well as the surface 
energies with the three corrections (green lines) are plotted for 
comparison. 

The prediction for the sticking velocity for dust ag gregates 



of 0.1 mg mass was approx. 10 ^ m s Mn the model of Giittlerl 



et al. pOlO ). The assumption behind that was that 95% of the 
kinetic energy in a collision is instantly dissipated due to plastic 
deformation (e ^ 0.2) and the remaining energy is balanced by 
the contact energy of the grains. However, a successful adapta- 
tion of the model of Thornton and Ning|| implies a sticking ve- 
locity in a more elastic regime, where sticking is dominated by 
adhesion rather than by plastic energy dissipation, i.e. the co- 
efficient of restitution is significantly larger than s - 0.2. This 
is clear from the mean values in Fig. [9] and also indicated by 
the green solid line. A sticking velocity of 10"^ ms"^ does not 
directly contradict our data but should instead be interpreted as 
a 50 % sticking probability. 

Recapitulated, the main simplifications were the following: 
(i) the average coefficient of restitution presented in Fig. [9] is 
also averaged over impact parameters (i.e., bjR ^ 0.5), which 
results in a systematic shift to higher values of s (cf. Fig[7|. Due 
to projection eff'ects, the unobserved component of the impact 
parameter makes this correction even greater. ( ii) The measure- 
ment of the coefficients of restitution is also constricted due to 
the projected view. It is reasonable to assume that collisions 
predominantly occur with major velocity components in the 
(observable) direction of the excitation (see Appendix B[ ). Af- 
ter the collision, the directions of the velocities are randomized 
due to the non-central collisions and irregular surfaces. This 
means that the velocity after the collision is more underesti- 



mated than the velocity before, which leads to a systematic un- 
derestimation for the coefficient of restitution. Qualitatively, 
this compensates the first shortcoming, while a deeper analy- 
sis is complicated. ( Hi) The eff'ective surface energy relies on 
the Hertz factor, for which we made an assumption of spheri- 
cal, monodisperse grains (logarithmic mean grain size) as this 
factor is impossible to describe analytically for irregular grains. 
(iv) The correction for porosity in Equation ^ assumes two 
porous surfaces, with only Hertzian deformation. Porous sur- 
faces are however rough and might interlink, which increases 
the eff'ective number of grains in contact. 

The sticking velocity given by IGiittler et al.| ( |2010| ) depends 
on the projectile mass as v oc ^-5/i« xhis is consistent with 
the model of [Thornton and Ning and also with our data for the 
larger dust aggregates. We therefore regard the current thresh- 
old of IGiittler et"an 

-5/18 



l^stick 



= (-) 



-5/18 



.-1 



m 
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(8) 



with m = 1.1 • 10"^^ kg as a good estimate for the velocity of 
a 50 % sticking probability. In Sect. |4.1| we showed that the 
thresholds for perfect sticking and perfect bouncing should be 
given by m = 1.9 • 10"^^ kg (sticking) and m = 2.1 • 10"^^ kg 
(bouncing). 

For particles with a given mass, this leads to a transition in 
velocity with a width of a factor of 325. The physical reason 
for this is not easy to deduce for these realistic particles, which 
feature a rough, irregular surface. The distribution of impact 
parameters is expected to be of major importance as it also has 
a large influence on the coefficient of restitution (cf. Fig. [7]). 
However, for a central collision we also expect a transition re- 
gion from perfect sticking to perfect bouncing if we consider 
a distribution of dust-aggregate sizes from 0.5 to 1.5 mm with 
an ellipsoidal shape. According to Eq. ^ smaller particles are 
sticking at higher velocities than larger ones. If those particles 
are oblate ellipsoids (aspect ratio 1:1.7:1.7, disk shaped) col- 
liding with their flat sides, the surface curvature at the contact 
site is enhanced by a factor of 1.7. For larger particles, we see 
from Eq. ^ that the sticking threshold velocity is generally 
expected to be smaller. If we furthermore consider two pro- 
late ellipsoids (aspect ratio 1:1:1.7, rod shaped) in the extreme 
case, which collide with the long axes aligned, the surface cur- 
vature is reduced by the factor 1.7. We can rewrite Eq. ^ as 
Vstick ^ R'^^^R^^^^, where the first term represents the mass de- 
pendence in Eq. (|4]), i.e., the inertia of the dust aggregates, and 
R' in the second term describes the influence of the surface cur- 
vature. The width of the sticking transition from the two efl'ects 
considered here is thus 

-3/2 / r,j^ ^2/3 

(9) 



vi 



IJRs 
^1/1.7 



where the indices s and 1 represent the small and large radii. 
This ratio can explain a width in velocity of a factor 5 for the 
transition from sticking to bouncing. This can be significantly 
larger if we also consider local changes in the radius of curva- 
ture (e.g., small bumps and dimples), which can be much more 
diverse. 
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Figure 1 1 : (color online) The image sequence shows a dust-aggregate dimer 
being hit by a single dust aggregate. First, the single dust aggregate collides 
with the upper dust aggregate of the dimer, which then rotates around the lower 
one. After this, the lower particle is hit as well and the dimer is accelerated to 
the right, but still intact. A movie of this collision can be found in the online 
version of this article. 



4.3. Contact strength of sticking dust aggregates 

As mentioned before, seven collisions led to sticking of the 
dust aggregates. The qualitative impression in the movie se- 
quences (see online material) is that the bonding between stick- 
ing particles appears rather weak. One of the dimers even 
collides with a glass wall (with unknown velocity), and the 
two particles separate again. However, one of the other dust- 
aggregate dimers collides with a single dust aggregate, and we 
have a good observation angle. This allows us to estimate the 
strength of the contact and obtain a quantitative estimate where 
the qualitative impression might be deceptive. We will first cal- 
culate a minimum shear strength of the connecting neck and 
then compare the kinetic energy in the collision with the adhe- 
sive energy in this neck. 

The image sequence of this interesting collision is presented 
in Fig. [TT] The dimer dust aggregate (particle 2) is nearly at 
rest and the single dust aggregate (particle 1) comes from the 
left colliding with the upper dust aggregate of the dimer (parti- 
cle 2a) at a collision velocity of 1.7 • 10"^ ms"^. The velocity 
vector is nearly perpendicular to the axis of the dimer dust ag- 
gregate. Particle 2a is hit, accelerated to the right, and then 
rotates around the common center with particle 2b. With a rel- 
ative velocity of 1.3 • 10"^ ms"^ particle 1 then collides with 
particle 2b of the dimer, causing the dimer to stop its rotation 
and move to the right from the point of collision. The two parti- 
cles of the dust-aggregate dimer are still attached to each other 
and can be followed until the very end of the experiment, 2.8 
seconds after r = s in Fig. [TT] during which the particles 
clearly move together. 

If the collision vector is perpendicular to the dimer axis, par- 
ticle 2a might be sheared off' if the momentum of the single 
dust aggregate is sufficiently large. We will calculate the shear 
strength of the connecting neck between particles 2a and 2b, for 
which this could happen, which is then a lower estimate for the 
actual shear strength of the dust-aggregate contact. If the neck 



has a cross sectional area A and particle 1 is pushing with a con- 
stant force F, the shear strength must be 5" > F/A to withstand. 
For the force we make an approximation to distribute the mo- 
mentum of particle 1 (mi • vi) equally over the impact time and 
assume a contact time of r = 5 ms as discussed by [Weidling 



etaL] ( [2Q09| ). With this we get 

miVi(l -h s) 



(10) 



For the contact area, we assume that the dimer was grown in a 
collision with elastic deformation and take the Hertzian contact 



area given by ( [Johnson[ |1985| ) as 



^Hertz = 



^ ^2 2\2/5 

15 m r^^'^ ^ 



(11) 



where nf and r* are the reduced mass and radius of particles 
2a and 2b, v = 4.3 • 10"^ ms"^ is the velocity of the collision 
from which the dimer had grown, and E = 8100 Pa is Young's 



modulus (see Sect. |4.2[ ). For the dust aggregates with radii of 
ri = 0.33 mm, = 0.25 mm, and r2b = 0.30 mm, this yields 
a contact area of Anertz = 9.3 • 10"^^ m^ and a shear strength 
of S > 640 Pa. This is a reasonable value and in agreement 
with compressive and tensile strengths in the range of 1000 Pa 
f Gu ttler et all [2009'). 

We can also compare the kinetic energy in the collision which 



is not dissipated due to plastic deformation (cf. Sect. 4.2) with 
the adhesive energy in the contact. For the kinetic energy we 
chose to take the reduced mass of only particles 1 and 2a as 
particle 2b appears to be nearly unaff'ected by the ffrst collision. 
The adhesive energy is N • £^break where N is the number of 
sticking contacts and £^break = 3.4 • 10"^^ J is the contact energy 
between two grains with radius a^ = 0.5 jim ( Poppe etal. 2000[ 
Blum and Wurm| |2000| ). With > 6:^^kin/^break we get > 
1700 contacts, which represents a contact area of 



Aadh - 



Nnal 



(12) 



With this, we get a contact area between the two dust aggregates 
Aadh > 3.8-10"^ m^, which is signiffcantly larger than calculated 
above. A more intuitive parameter is the contact radius which 
is ^^Hertz = 17 yum and ^adh = 35 yum, respectively. 

We have to note that the energy consideration has the prob- 
lem that it assumes that all the kinetic energy left after plas- 
tic deformation goes into breaking the contacts. For the shear 
strength considerations, we showed that the velocity to break a 
dimer is proportional to the shear strength, which may well be 
a factor of two higher than the minimum calculated here. This 
yields a breakup velocity around 0.3 m s"^ which is higher than 
the highest sticking velocity found in our experiments (see, e.g.. 
Fig. [8|. Thus, if we ffnd conditions that favor sticking in a qui- 
escent region in a protoplanetary disk, sticking will continue 
and further collisions at the same velocity will not break con- 
tacts again. 

4.4. Consequences for the protoplanetary dust growth 

The consequences of our experiments for the growth of pro- 
toplanetary dust aggregates are not easy to predict. The model 
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of |Guttler et al.| ( [20TQ| ) is highly non-Hnear and small changes 
can have significant impact. However, it is possible to give an 
indication about the mass of the dust aggregates (relative to the 
mass computed by |Zsom et aL] |201Q[ |2011| ) at the time when 
the dust aggregates cross the boundary between sticking and 
bouncing for the first time. 



In Sect. 4.2 we showed that the sticking threshold velocity 
scales with mass according to Eq. [8] To determine the ratio of 
two masses for two diff'erent scaling factors of the threshold line 
mold ( |Guttler et al. 2010, ) and mnew (this work), we can write 



^old 



18/5 ~ 
^new 1 ^new 



Void 



^old 



(13) 



where the first term on the rhs. denotes the ratio of typical 
collision velocities for the dust aggregates with the respective 
indices. If the relative velocity for the dust aggregates in the 
protoplanetary disk was not dependent on their size, this term 
would be unity and we just get the ratio of the two scaling fac- 
tors. However, we know that the relative velocities are increas- 
ing with the dust aggregates' masses as v oc m^^^ under the 
conditions considered here ( [Weidenschilling and Cuzzi 1993 ), 
and thus we get 



•^new 
Void 



^old 



1/3 



Using this for Eq. 13 we arrive at 



Wold 



Wold 



5/11 



(14) 



(15) 



which yields a significant factor of 250. 

An exact mass for the largest dust aggregates that can di- 
rectly grow shall not be given here, as this number could easily 
be misinterpreted. It must be stated that this is the enhanced 
growth only for one certain time in the evolution of protoplan- 
etary dust. After this, new non-linear eff'ects could show up in 
the model of [Giittler et~aL] ( |2QT0| ) and |Zsom et al.| ( |2QTQ| so 
that the overall growth might be enhanced in the same way but 
could as well stall at the same mass as before. An interesting 



aspect is, however, that the so-called 'bouncing barrier' (Zsom 
et aH |2010) is becoming smaller: growth is enhanced and also 



the fragmentation threshold was recently revised (Beitzetal 



2011 



and found to be at smaller velocities for large (centime- 
ter) dust aggregates. Thus, it might be easier to trigger fragmen- 
tation, widen the dust-aggregate size distribution, and by this 
open a new path for the further protoplanetary dust- aggregate 
growth. 



5. Conclusions 

With a novel experimental setup, we conducted micrograv- 
ity experiments with a multitude of porous dust aggregates of 
diff'erent sizes and analyzed free collisions between them. We 
found a total of 125 collisions with relative velocities between 
2 mm s"^ and 1.7 m s"\ where 22 collisions happened between 



1.0 to 2.0 mm-sized dust aggregates that all resulted in bounc- 
ing and 103 collisions between 0.5 to 1.5 mm-sized dust aggre- 
gates (see Fig. [5]). Of those 103 collisions we observed stick- 
ing in seven cases, bouncing in 95 collisions and once one of 
the collision partners fragmented. This was the first time that 
sticking of millimeter- sized dust aggregates was observed un- 
der controlled experimental conditions. 

The sticking collisions found in the experiment scatter over 
a wide range of velocities, all of which are higher than pre- 



dicted by the model of Guttler et al.| ( [20TQ| . Therefore, we sug- 
gest to implement a transition zone from sticking to bouncing 
collisions with a velocity-dependent probability for the colli- 
sion partners to stick to each other. The threshold line given 
by IGiittler et al. is consistent with a 50 % sticking proba- 
bility in our experiments, while the velocity thresholds for a 
100% and 0% sticking probability are given by our Equation 
^ with m = 1.9 • 10"^^ kg (100% sticking probability) and 
m = 2.1 • 10"^^ kg (0% sticking probability). 

To determine a general description for the fastest velocity 
that always leads to sticking, we adapted the collision model by 
[Thornton and Ning| ( |1998| ) in Sect. |4.2[ We determined Young's 
modulus of our dust aggregates to be £^ = 8100 Pa by deriving 
a yield velocity of Vyieid ~ 9 • 10"^ ms"^ from the decline of 
the coefficient of restitution with increasing collision velocity. 
By accounting for porosity of the dust aggregates and the Hertz 
factor we received a term for the effective surface energy of our 
dust aggregates, which we use to formulate an analytic term for 
the sticking velocity as given by Equation (|4]). 

Furthermore, we could show that the contact between stick- 
ing dust aggregates of millimeter size is sturdy enough to let 
the dimer survive collisions with single dust aggregates if the 
velocity is low enough to provide a chance for sticking. 
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Appendix A. Gas drag on the dust aggregates 

Particles moving in vacuum under microgravity conditions 
are expected to do so on linear trajectories with constant veloc- 
ity. However, disturbances like static charging of the particles 
or drag of the residual gas may alter this behavior. While we 
can rule out the former for those 1 mm dust aggregates undergo- 
ing collisions by having checked their trajectories for parabolic 
components we have to address the latter more thoroughly, as 
gas drag slows down the particles exponentially even without 
collisions. A measure for the importance of gas drag is the 
stopping time Tf , after which a particle's velocity has decreased 
to 1 le of its starting value. With a gas pressure of p = 10 Pa the 
mean free path of the gas is 1 mm, which is comparable 
to the radius r of our particles. Therefore, we calculate the fric- 



tion time both for the Stokes and the Epstein regime after Blum 
( ,2006, their Sect. 5.1). They give the stopping time as 



m (Epstein) 



with pb = 2600 kgm (Blum et al. 
1.2 kgm"^ • {pIlQ^ Pa), 
Vg = 655 ms 



2006 



Table 1), Pg = 
17.2 yuPas and 



= 0.347, //g = 
^ being the bulk density of the dust particles, 
the gas density, the filling factor of the particles, the viscosity 
and the thermal velocity of the gas, respectively. With par- 
ticle radii between = 0.25 mm and rmax = 0.75 mm, 
this corresponds to stopping times between 0.7 and 9.3 s, de- 
pending on the particle size and the friction regime. Although 
both times are of the order of the experiment duration, most 
trajectories are rather short due to inter-particle and particle- 
wall collisions and the exponential function can be approx- 
imated linearly over short sections (here 60 data points, i.e. 
120 ms, were used at most). The frictional velocity decay 
is thus between Av/v = 1 - exp(-0.12 s/0.7 s) = 0.16 and 
Av/v = 1 - exp(-0.12 s/9.3 s) = 0.01 for tracks that have the 
maximum length. Most of the tracks were shorter than that, 
reducing the frictional decay even more. 

Appendix B. Projection effects 

In order to investigate the error we introduce by using just 
one projection and not the three-dimensional data we compared 



14 




velocity [m/s] 

Figure B.12: The velocities in the x- and z-direction were derived directly from 
the fits. The velocities in the shaking direction are higher than perpendicular to 
it. The median velocities are Vx = 2.9 mm s~^and = 4.6 mm s~^ . 



the components of the particle velocities. The z-component of 
the velocity, pointing parallel to the shaking movement, is the 
most prominent term. While in the other directions the parti- 
cles reach kinetic equilibrium, the shaking induces additional 
energy in the z-direction. Figure B.12 shows the fraction of 
the particles with a certain velocity or smaller throughout the 
whole experiment. It can be seen that the particles on average 
are faster in the z- than in the x-direction. The median veloc- 
ity in the x-direction is 2.9 • 10"^ ms"^ while it is 4.6 • 10"^ 
ms"^ in the z-direction. We do not expect the velocities in 
the }7-direction to differ from those in the x-direction, so that 
the average difference between projected and total velocity is 
V4.62 2 • 2.92/ V4.62 -h 2.92 - 1 - 0.13. This would shift the 
results a little, but not change the overall picture. 
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